Radiation of the atomic chains exited by a channeling particle 



V. Epp and M.A. Sosedova 
Tomsk State Pedagogical University 

Abstract 

(N 

Basic properties of radiation of the atomic chains excited by a channehng particle are considered. 
Using a very simple two-dimensional model of a crystal lattice we have shown that the main part of this 
radiation is generated on the frequency of oscillations of a channeling particle between the crystal planes, 
shifted by the Doppler effect. Angular distribution of the radiation of the chain of oscillating atoms is 
sharply peaked in the direction of the velocity of channeling particle because of coherence of the fields, 
produced by individual atoms. 

^ 1 Introduction 

O Physics of channeling of the accelerated particles in crystal, is a quickly developing field of science [T]. 

^ The effect of channeling has served as the base of development of new experimental methods of research 

on the crystal structure. By use of the channeling effect it is possible to study thermal oscillations and 
^ displacement of atoms in a lattice, distribution of electronic density in interatomic space of crystals and 

• ^ to measure exact orientation of the crystal planes. Channeling of light particles is used for generation of 

intensive monochromatic X-ray radiation. Using experiments and the theories which describe the orientation 
effects in crystals, a new sources of X-ray and the gamma radiations consisting of accelerators and precisely 
Qh oriented crystals are created. 

In the course of interaction between the channeling particle and a crystal the electromagnetic radiation 
^-H is generated. In the framework of classical electrodynamics one can consider this radiation as radiation from 

^ different sources - radiation of the channeling particle oscillating between the atomic planes or around an 

crystal axis, radiation of the electronic gas excited by a channeling particle (weik fields of electrons) and 
the radiation which is produced by the excited atoms of a crystal lattice. Radiation of a channeling particle 
is investigated in details (see for example [21 [3]), the experimental research of radiation of weik fields can 
be found in f?|. Radiation of the atoms excited by a channeling particle is has not been studied so far. 
In this article we consider part of this radiation caused by oscillations of the atoms embedded at specific 
sites in the lattice of a crystal, leaving aside the radiation caused by quantum transitions of electrons from 
excited to principal states. We suppose that the conditions of implementation of classical electrodynamics 
are fulfilled. Namely, the energy of excitation of atom is much more than distance between the energy levels; 
we consider the part of a radiation spectrum of single atom where the energy of photons is much less than 
energy of excitation of atom. Specific feature of considered radiation is that the phases of oscillations of 
atoms correlate between themselves as oscillations are excited by the same channeling particle. Therefore 
radiation should be considered as coherent. 
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2 Dynamics of atoms in a crystal lattice 

Let us consider very simplified model of a crystal lattice which allows, nevertheless, to find out the basic 
properties of radiation of the excited atomic chains. We will consider a crystal as a two-dimensional lattice 
of right-angled structure. Let the relativistic, positively charged particle having a charge qi enter in a crystal 
between two atomic chains under small angle to the crystal axes. We consider the energy of the particle 
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much above the energy loss on radiation and interaction with the crystal. Then the velocity of the charged 
particle in the crystal remains to be constant and equal to the initial velocity V. 

Now we are interested in oscillations of the atoms at the lattice points which is caused by collision with 
the channeling particle, and in the radiation related to these oscillations. In spite of the fact that the atom, 
as a whole, is neutral, only part of its electrons are bound to the atom. The peripheral electrons belong 
to the whole crystal. Colliding with an atom, a channeling particle interacts with the screened potential of 
the atom nucleus, disturbing it from the equilibrium. Oscillations of a nucleus with interior electrons causes 
time-dependent polarization of atom and hence electromagnetic radiation. At calculation of interaction of a 
particle with individual atom we assume that the channeling particle moves parallel to atomic chains with 
a constant velocity. Rather slow oscillations of a particle between atomic planes we will take into account 
by the distance between a particle and nucleus of atoms, which is slowly changing along a trajectory of 
channeling particle. The model of a crystal lattice considered here has been used in [5 for calculation of 
polarization of atomic chains induced by the channeling particle. 

Let the axis OX of a coordinate system be parallel to the atomic chains and lying in the middle between 
the neighboring atomic chains of a crystal (fig. [T]) . We denote the distance between the next atoms lying in 
a chain by b and the distance between the atomic chains by 2D. 
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Figure 1: Model of a two-dimensional crystal lattice. 



Then the atom with number i of the first chain has coordinates {xi,D), and the opposite atom of the 
second chain has coordinates {xi, —D). We choose the origin of coordinates system so that xq — 0. 

The field of a relativistic channeling particle in the laboratory reference system is concentrated basically 
in a plane which is perpendicular to a direction of motion [6]. 

Hence, the particle interacts effectively with the atom during a short interval of time when the atom is 
in the vicinity of this plane. One can expect that this interaction will be reduced to the momentum transfer 
from the particle to atom in a direction orthogonal to the velocity of channeling particle. In order to prove 
it, we find x and y components of the momentum transmitted to the atom of a lattice. The electric field of 
a charged particle moving with constant velocity is set by the formula [6J: 
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where ^ is a vector pointing from the charged particle to the atom, /? = V/c, c is the velocity of light, '0 is an 
angle between the velocity of the moving particle and vector R. Vector n is a unit vector of direction defined 
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by n = R/R. The angle ip depends on parameter a which is the shortest distance between the trajectory of 
channehng particle and a nucleus of atom. 

a a 

simp = — = (2) 
R + VH^ 

Time t is equal to zero when the particle is at the minimal distance from a nucleus. The atom is acted on 
by the field of moving positively charged particle with the Coulomb's force 

F = Eq2. (3) 

where q2 is an effective charge of a nucleus of atom which is screened by the interior electrons. It is obvious 
that atoms of a crystal lattice are neutral, but a channeling particle approaching the lattice point of a 
crystal appears to be inside the atom and interacts therefore with the nucleus partially screened by the 
interior electrons. Consequently, an effective charge (72 is less than the charge of a nucleus. 

One can see from the formula ([T]) that the field of an ultrarelativistic particle (/3 '-^ 1) has important 
quantity in a vicinity smip ~ 1 or more exactly cosV' ~ 7~^, where 7"^ = -^/l — is the relativistic factor. 
Then, it follows from ([2]) that interaction of a particle with a nucleus occurs during time: 

a a 

Integrating coordinate components of formula ^ over time i G [— cxd, 00], we find components of momentum 
Px and Py which is transmitted to atom of chain along the axes OX and OY respectively: 
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p. = f|^yr3^, (5) 



It follows from the ratio 



p. 



P. = ^. (6) 



= 1, 



that the momentum transferred to the atom in Y direction is much greater then that in X direction. Hence, 
in case of an ultrarelativistic particle motion of the atom along the axis OX can be neglected and only 
oscillations along the axis OY can be considered. The passing particle excites harmonic oscillations of the 
atom at frequency uj which is defined by crystal properties. For the moment we consider oscillations of atom 
as free ones, attenuation of these oscillations we shall take into account later. The law of motion of the i-th 
atom of the first chain can be written as: 

T .T ^ ib\ lb 
yl = A'^sinu[t--J+D, t>-. (7) 

Accordingly, the law of motion for the atom of the second chain will be written as: 

.2 _ a2- (■ ib\ „ ib 



yt^-Afsmujl^t--j-D, t>-, (8) 

where {j = I, 2) is a constant defined by initial conditions. The origin of coordinates is accepted so that 
t = when the channeling particle hits the atom number i = 0. The amplitude of oscillations is defined by 
an impulse transmitted to atom. It is obvious that the laws of motion Q and ([s]) should be completed by 
equation 

i ib 

yi^o, t<y- 

Now, in order to calculate the amplitude of oscillations we suppose that the particle transmits part of its 
momentum to the atom in no time. Such assumption is valid, if effective time of interaction of a particle with 
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an atom is much less than period of its own oscillations. From the formula Q follows that this condition 
takes form 

uja 

— <7, 
c 

which obviously is fulfilled for a nonrelativistic atom as ioa is of the same order of magnitude as the velocity 
of atom oscillating with amplitude a. The amplitude of oscillations we find from initial condition that the 
amplitude of momentum of the atom raujA'^ is equal to momentum (j6| accepted from the particle. As a 
result we have: 

mLuV a\ 

where is the shortest distance between the particle and i-th atom of chain number j. As the velocity of 
atom is the nonrelativistic, one can use dipole approximation at calculation of radiation of an oscillating 
charge. We will combine two atoms which are in different chains opposite each other and are simultaneously 
excited by a passing particle into one dipole. The dipole moment of obtained system is equal to Pi = 
(2xi(7i, pyi, 0), where 

Py,{t) = qM] - Al) sinu (t-'^y (10) 

Let us find dependence of parameter on coordinate x. Motion of a channeling particle in case of small am- 
plitude of oscillations (it means the amplitude of oscillations in Y axis direction is much less than interatomic 
distance 2D) to a good approximation can be described by the equations: 

yc = K cos{nt + (f>o), Xc = Vt, 

where j/c and Xc are particle coordinates, K and 00 is amplitude and initial phase. The initial phase can be 
putted equal to zero by corresponding choice of the origin of X axi^ Then the equation of trajectory of a 
channeling particle takes the form: 

riXf^ 

Vc = Kcos—. 

Accordingly, parameters aj for the atoms having coordinates Xi = ih, are equal to 

a\=D~Kcos^, a^ = D + Kcos^. (11) 



Substituting these parameters in formulas (|9| and (10 1, we obtain amplitude of oscillations of i-th dipole 

A. = ,,(4-AJ = ^^^^cos— , (12) 
And the net result for the dipole moment generated by two atoms of the neighbor chains will be written as: 

Pyi{xi,t) = AiSiiiui {t-^^ , Pj:i{xi) = 2x,qi. (13) 



In principle, the form of expression ( 13 1 was obvious from the beginning. All previous reasoning has been 



directed to calculation of amplitude of oscillations (12). At this point we can correct the primitive model of 



a crystal lattice to some extend setting the amplitude Ai from some empirical data. For example, we have 
not set the quantity of screened charge of nucleus q2, have not considered thermal oscillations of atom etc. 



Some of such corrections can be made by inserting of corresponding coefficients in formula (12 1. 

The most essential improvement we are going to make at this step - is to take into account the attenuation 
of atom oscillations. The main cause of attenuations is the transmission of the energy of oscillations to 
surrounding atoms and to a lattice as a whole. In terms of quantum mechanics this process is described by 
radiation of phonons. By assuming that the energy of phonon is much less than the energy of oscillations 



^We have already defined the coordinate origin just after Eq. (jsjl, bounding it to some specific atom. Now we can shift the 
origin by integer number of distance b. Note that b <g; V/Q. 
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we can consider attenuation of oscillations by putting an exponential factor into amplitude of oscillations. 
With account of attenuation the equation for the dipole moment ( 13 1 becomes: 



A,e-"(*-^)sina; (^i- ^) , t> 



Py^i^^,t) = i ' ^ VI// V (14) 

. ' - y 

where a is the attenuation coefficient. We will ignore the constant component of the dipole moment Pxi as 
the radiation field is proportional to the second derivative on time from the dipole moment. 



3 Radiation of the excited chains of atoms 



Electric field of radiation of the dipole moment is defined by the formula [5] : 

1 



(15) 



where the unit vector n is defined by equality n — f/r, r is the vector pointing from the dipole to the point 
for which the vector Ei is written down, r = \r\, Pi is to be evaluated at the retarded time t' 



t' = t 



Further we assume that the distance between the crystal and the observer is much greater than sizes of the 
crystal. Then the vector n and distance r in Eq. (15 1 are constant and do not depend on x. Dependence 



on X should be considered only in the phase of an electromagnetic wave, namely in t' . We denote the vector 
from the coordinates origin to the point of observation by and radius- vector of the dipole moment by Vp 
as shown in Fig. [2] As i? 3> r^, it is possible to write approximately: 



Y 




where n — R/ R. Accordingly, for the retarded time t' we have: 



t' ^t- 



R — xuj 



(17) 



Let us find the components of vector E{t) in a spherical system of coordinates shown in Fig. [sj Direction of 
the unit vector n is defined by the polar angle 6 and azimuthal angle tp. The Cartesian coordinates of the 
unit vectors of spherical coordinate system can be written as: 




Figure 3: Spherical system of coordinates. 



Then it follows from Eq. (15) that: 



(sin 9 sin (p, cos 6, sin 9 cos tp) , 
(cos 9 sin(/3, — sin 9, cos 9 cos (p), 
(cos ip, 0, — sin ip). 



E^{R,t) = ^py,{t') sin t 



AiC [(a^ — w^) sinwT — 2q!WC0swt] 



Taking twice derivative from Eq. (14) with respect to time, we obtain 

Pyi (-^z ; t ) 

where 



0, t' < 



t' > 



v 



(18) 
(19) 
(20) 

(21) 



(22) 



r = t'-5=t-^-5(l-/3n.), /3=^. 
V c V c 

Let the atomic chain consist of large, but finite number of atoms N. The net field of radiation is the discrete 
sum of fields of individual atoms of the crystal 



N 



^(i?,i) = ^^,(i?,i). 



(23) 



As the distance between the atoms is much less than the wavelength of radiation, it is possible to turn to 
continuous distribution of the dipole moment along the X-axis, and replace the sum by integral. There is 
one dipole per interval of length b in atomic chain, hence, the quantity e{x, t) = Eie{R, t)/b has the meaning 
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of the field generated by unit lengtli of the atomic chain. Substituting into this expression the field of Eq. 

i: 

(24) 
(25) 



(21 1 and taking into account Eq. (22 1, we obtain e(x, t) for a time interval where it is not equal to zero: 

e{x, t) = Aoe~°''^ cos [{a^ — iJ^) sinwr — 2au) cos wr] 



here 



2qiqlKsine 

Aq — „ ■ „ „ , T = t 



R X 



(l-/3n.). 



mVD^ujbRc^' " c y 

Trigonometric functions in the square brackets can be combined, using the phase ^p•. 



sin ^uj 



2 ' 



COS i/juj 



As a result we get: 



fix 



e{x, t) = Api{a + Lo )e "'^ cos sinti;(r — ip) 



Integrating this expression over x, we find the field generated by pair of atomic chains 



E{t)6 



e{x, t) dx. 



(26) 
(27) 

(28) 



The lower limit xi of integral is coordinate of the crystal origin. The upper limit X2 is a time function. 

While the channeling the particle moves in a crystal, the upper limit coincides with coordinate of the 
particle according to Eq. (22). After the particle leaves the crystal, radiation is generated by all the chain 
of atoms, and X2 coincides with coordinate of the end of chain 0:20: 



X2{t) 



[l-pn^) c V 



R X 

X20-, t> — + —(1 - I3n^). 
c V 



(29) 



Integration in Eq. ( 28 ) gives 

Ee{t) = ^Ao{a^ + uj^) e 



where 



tn = t- 



R 



■a{to-kx) ( asin(a;ti + fca;A_) - A_ cos(a;ti + fca;A„) 

a2 + 

asin(a;f;i — fc.TA+) + A+ cos{ujti + fc.TA+) 
a2 + Ai 



R 1 

ti=t V, fc=-(l-/3n^), A±=l]'±w. 

c V 



(30) 



(31) 



We see from this expression that radiation is generated on two frequencies, namely on frequency of atoms 
oscillations lo and the frequency of oscillations of channeling particle which is shifted by the Doppler effect: 



Eq. (30) can be written as: 



Eg{t) = e-"*° [F(a;2,t)e"'="^ - i^(:Ei, t)e"*^"i] , (32) 
where F{x,t) is periodic function of time. Considering exponential factor in last formula, we see that while 

R X , ^ , 

c V 
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the amplitude of the electric field of radiation is growing accordingly to 

and approaching a constant level if the crystal is long enough. It means that the time it takes for the particle 
to cross the crystal is much greater then the lifetime of oscillating atoms 1/a. We shall refer such a crystal 
as a "long" one. After a particle leaves the crystal, the amplitude of field exponentially decreases. Let us 
next neglect the boundary effects and consider only radiation with constant amplitude of the field. In other 



words, we consider only first term in Eq. (32) where X2 — t^/k. This is the main term in case of a "long 
crystal" . Then the formula for a field of radiation becomes 



Ee{t) 



I 

2k 

cos ri't 



2) <^ sinf^'t 



a cos Loip — IS. ^ sin Loip a cos wi/' + A+ sin 



— asinw?/; — A_ coscj?/' a sin cj^/; — A_(- cos w0 
a2 + A2 a2_^A2 



(33) 



As we see in this case the radiation is generated only on frequency fi'. Let us find angular distribution of 
intensity of radiation. Intensity of radiation dl in an element of space angle do is defined by the formula 



dl 
do 



An 



Averaging expression ( 33 ) over the period of radiation and substituting in ( 34 1 we obtain 

dl _ qjqjK^ sin" 9 AQ'^a^ + ja^ + uj^)^ 



(34) 



(35) 



From last formula follows that dependence of intensity of radiation on frequencies f2' and w has resonant 
nature. At a reasonably small a the resonance occurs at fi' = a; or 



(36) 



If the channeling particle is ultrarelativistic one, sufhcient wide range of values of the relation w/fi' satisfies 
the resonance condition (36): 

1 ^ „ 9 



4 Conclusion 



In this work the basic properties of radiation of the atomic chains excited by a channeling particle are 
considered. On a simplified model of a crystal lattice it is shown that this radiation is generated on two 
frequencies, namely on the frequency of atom oscillations, and on the frequency of oscillations of a channeling 
particle, shifted by Doppler effect. In case of a reasonably long crystal the radiation with frequency of 
oscillations of a channeling particle prevails. Occurrence of the factor (I — ^n^) in denominator of the 
formula for angular distribution of intensity of radiation ( 35 1 results in relativistic effect - the basic part of 



radiation is generated into a cone around the direction of velocity of a relativistic channeling particle. Hence, 
the frequency of radiation and angular distribution are similar to corresponding characteristics of radiation 
of a channeling particle. This can rise a problem of experimental separation of one form of radiation from 
another one and a problem of identification of considered radiation. 

The interesting effect showing that a set of not relativistic oscillators at rest or oscillating atoms produces 
radiation with the properties which are specific for radiation of the relativistic charged particle, is obviously 
a result of the fact that oscillations of atoms are coherent. Because phases of their oscillation are tuned by 
the channeling particle. The interference of fields of radiation of individual atoms gives the properties of the 
net radiation, described above. 
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